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Abstract. We present results of anisotropic thermal expansion and low temperature 
magnetostriction measurements on YNi2B2C single crystals grown by high temperature 
flux and floating zone techniques. Quantum oscillations of magnetostriction were 
observed at low temperatures for H\\c starting at fields significantly below Hc2 
{H < 0.7 Hc2)- Large irreversible, longitudinal magnetostriction was seen in both, in- 
plane and along the c-axis, directions of the applied magnetic field in the intermediate 
superconducting state. Anisotropic uniaxial pressure dependencies of Tc were evaluated 
using results of zero field, thermal expansion measurements. 
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1. Introduction 

The members of the i?Ni2B2C {R = Gd-Lu, Y) series of compounds serve as model 
systems for studies of a number of phenomena: co-existence of local moment magnetism 
and superconductivity, non-locality and flux line lattice transitions, heavy fermion 
physics and complex metamagnetism [U El El El- Although many sophisticated 
experiments were performed on i?Ni2B2C borocarbides, data on thermal expansion 
(TE) and magnetostriction (MS) for several members of the series exist only for limited 
temperature ranges, or only on polycrystalline samples, or are not available at all. For 
one of the less complex members of the series, non-magnetic, superconducting YNi2B2C, 
for example, only a rather limited set of results has been disseminated [3 EI. In 
this publication we report anisotropic (c-axis and a6-plane) TE and longitudinal, low 
temperature, MS measurements on YNi2B2C single crystals grown by two different 
techniques and well-characterized by other methods. The objective for this work is 
manifold: to acquire data that can serve as a baseline in studies of TE in more complex, 
magnetic or strongly correlated, borocarbides; to probe the irreversible properties 
in intermediate superconducting state; and finally, to evaluate the uniaxial pressure 
derivatives of YNi2B2C from the anomalies in TE at Tc (using the the thermodynamic 
Ehrenfest relation), an approach proven viable e.g. in HTSC fHl, 1^, and compare the 
results with the implications of the analysis in Ref. 

2. Experimental methods 

Two YNi2B2C single crystals were used in this work. One of them (sample A throughout 
the rest of the text) was grown by a Ni2B high temperature flux method (see Ref. 
m El E] for more details). This method yields plate-like crystals with the c-axis 
perpendicular to the plates. The crystal was shaped into a nearly rectangular bar, with 
two pairs of parallel surfaces, so that TE and MS were measured along [100] (L = 2.37 
mm) and [001] (L = 0.92 mm) directions. Another crystal (sample B) was grown using 
vertical zone melting method with a commercial 4-mirrors image furnace (model FZ-T- 
4000-H-VI-VPM-PC, Crystal Systems Corp., Japan). We used drop cast polycrystalline 
rods (~ 6 mm diameter, ~ 60 mm long) as feeding rod and crystal support. High purity 
argon at atmospheric pressure was used as a protective atmosphere inside the quartz 
working tube. Feeding and bottom shafts were rotated in opposite directions to insure 
an effective stirring of the molten zone, typically at +10 and -35 rpm. During the 
growth, the growth speed was decreased from 10 to 2 mm/h to favor the formation 
of large, single grains. The phase purity of the single crystals was checked by Debye- 
Sherrer powder diffraction pattern and compared to original feeding rod spectrum. No 
detectable extra reflections due to impurity phases or changes of cell parameters were 
observed. An as-grown bar was oriented with a Laue camera and cut into a close to 
prismoidal shape with two pairs of parallel surfaces, so that TE and MS were measured 
along [110] {L = 3.03 mm) and [001] (L = 1.82 mm) directions. After shaping, both 
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samples were annealed in dynamic vacuum (10"^ — 10~^ Torr) at 950° C ^S]- It should 
be mentioned that for reasons unrelated to the objectives of this work, the sample A 
was grown with isotopically pure ^''B and the sample B with ^^B. 

Magnetization measurements were performed using a commercial MPMS-5 
(Quantum Design, Inc.) SQUID magnetometer; ac resistance (/ = 16 Hz, / = 3 mA) 
and heat capacity were measured with ACT and heat capacity options of a PPMS-9 
instrument (Quantum Design, Inc.). 

Thermal expansion and magnetostriction were measured using a capacitive 
dilatometer constructed of OFHC copper; a detailed description of the dilatometer 
will appear elsewhere Plj. The capacitance was measured with an Andeen-Hagerling 
2500A capacitance bridge (the bridge resolution of 10"'' pF corresponds to a sample 
dilation measurement limit of about 0.3 pm when the dilatometer is operating near 
20 pF). The dilatometer was mounted in a PPMS-14 instrument (Quantum Design 
Inc.) and was operated over a temperature range of 1.8 to 305 K and in magnetic 
fields up to 140 kOe (either in vacuum or with a small amount of helium exchange 
gas to minimize thermal gradients). The temperature was measured with a calibrated 
Cernox-1030 thermometer (Lakeshore Cryotronics) mounted on the dilatometer. The 
field dependence of the dilatometer was less than ±3 A over ±140 kOe and will be 
ignored. The temperature dependence of the dilatometer (or "cell effect") is removed 
using published values of the thermal expansion of copper ^H] • Data were acquired with 
the temperature increasing at a rate of about 0.4 K/min or with the field changing at 
a rate of about 6 kOe/min. The absolute accuracy of the dilatometer was checked with 
measurements on a 4 mm sample of pure aluminum that were compared to published 
values ^H]- The average deviation between our aluminum TE measurements and those 
in the literature over our full temperature range is 7.0 x 10^^ K^^. The average fractional 
deviation of our aluminum measurements from those in the literature is about 1% above 
40 K but becomes larger at low temperatures where the thermal expansion of aluminum 
is sensitive to both sample purity and preparation ^^l- The maximum deviation of 
2.8 X 10~^ occurred at 300 K and corresponds to a fractional deviation of 1.2% 
from published values. We estimate an experimental uncertainty of 10~* in our 
thermal expansion measurements from the mean deviation to a fit of our aluminum 
data near 20 K, this corresponds to a dilation uncertainty less than 0.5 A/K. 

3. Results and Discussion 

3.1. Thermal expansion 

Superconducting transitions, as measured by DC magnetization, zero field resistivity 
and heat capacity, are shown in Fig. ^ In all three measurements the transitions 
are sharp, being sharper for sample A. The Tc values as defined by p = 0, onset of 
diamagnetism in M{T) measurements and mid-point of the heat capacity rise (or by 
balancing of the normal state and superconducting state entropies) are ~ 15.9 K and 
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~ 15.1 K for samples A and B respectively. This difference in between samples A 
and B is partially due to the boron isotope effect [121 Ej and partially due to different 
scattering in two samples as can be see from the difference in the residual resistivity 
ratio, RRR = p(300K)/p(17K) ^ 63 and 16 for A and B respectively [12] . It should be 
noted that sample B had a small, rather broad, bump in resistivity near 30 K, feature 
that was not observed in the sample A or other fiux-grown YNi2B2C samples. 

Temperature-dependent linear and volume (defined here as AV/Vq = AL^/ Lco + 
2ALab/ Labo) dilations (relative to the values at 1.9 K) for two samples are shown in Fig. 
121 The two crystals behave very similarly with the slight discrepancy most probably 
reflecting the accuracy of the measurements of the room temperature dimensions and 
their changes and possible, slight imperfections in the shape of the samples. In-plane 
linear and volume thermal dilations are positive in the temperature range studied. The 
c-axis thermal dilation is negative at low temperatures, it changes its sign at about 200 
K. 

Temperature-dependent linear and volume thermal expansion coefficients measured 
on samples A and B are shown in Fig. jSl cxdT) data for both samples are very similar, 
as well as the aabiT) data above ~ 100 K. Further studies will be required to understand 
whether the broad maximum in a^ioo] {T) of the sample A at about 80 K is a real feature 
albeit it was reproducible in two consecutive measurements:, one in high vacuum mode 
of the PPMS, another with low pressure He exchange gas in the sample chamber. Our 
data are compared with the published data [HI H] obtained by powder and single crystal 
X-ray diffraction at different temperatures. To be consistent with the Fig. 3 and results 
for (3 in the Table 3 of Ref. [7|, the data cited in ref. [7| for YNi2B2C should be 
divided by 10. Such "corrected" Oc data together with the aa and [3 data from the Table 
3, Ref. [7] are included in Fig. (21 Our data are consistent with the published values 

[ilZj. 

The temperature dependence of the linear thermal expansion coefficients close to Tc 
is presented in Fig. jU Jumps in aiiT) are seen unambiguously, with the opposite sign 
of the jumps for aab and ac and very similar, within the errors of the measurements, 
values for both samples. Using results from the heat capacity (Fig. and thermal 
expansion coefficients (Fig. ^ measurements, the uniaxial pressure derivatives of the 
superconducting transition temperature, Tc, can be calculated using the thermodynamic 
Ehrenfest relation: 



^ = Vmoi ■ Aai{Tc 
dpi 



ACpjT, 



where Vmoi is a molar volume of the material (for YNi2B2C, Vmoi = 39.5 cm^) and 
Aaj(Tc) and ACp{Tc) are the jumps in the i-th thermal expansion coefficient and specific 
heat at T^. Estimates of the uniaxial pressure derivatives for YNi2B2C based on the 
aforementioned measurements are summarized in Table ^ Within the 10-20% error 
bars in Aoj (see Fig. there is no significant difference in estimated uniaxial pressure 
derivatives for two samples. In the table dT^/dP* = 2-dTc/ dpab + dT^/ dp^.. The dT^/dP* 
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defined in such way lacks the contribution from the off-diagonal dT^/ dpij terms that play 
a role in the experimentally measured dT^/ dP under hydrostatic pressure. Since the off- 
diagonal terms are usually significantly smaller then the diagonal ones, we can still 
compare the last column of the Tabled with the experimentally measured hydrostatic 
pressure derivatives. The experimental data on dT^/dP of YNi2B2C were presented in 
a number of publications, some of the reported values are listed below (in K/kbar): 
-0.58 ■ 10-2 [TH], -0.9 ■ 10-2 [19J, 0.32 ■ 10"^ (in P ^ limit) [20], and -0.9 ■ 10"^ 
|21j . All but one of these results are consistent with our estimates. We are unaware of 
any direct uniaxial pressure measurements on YNi2B2C other than very short conference 
proceedings publication j22] , that contains a statement that for this material is almost 
independent on (uniaxial) pressure, again in agreement with our estimates of dT^/dpi 
(Table IT}. 

The uniaxial pressure derivatives obtained from the Ehrenfest relations in this work 
differ in value and, more importantly, have the signs opposite to the ones inferred in 
Ref. JO]- The procedure used there was indirect, based on a number of assumptions 
(the most important, but not universally correct, being equivalence of physical and 
chemical pressure), had to rely (due to the nature of the doping used) on refining of 
the structural data of two-phase samples, and therefore is open for discussions and 
criticism. On the other hand, it is noteworthy that additionally our current results 
(Table apparently contradict the implications of the band-structure calculations 

that superconductivity in RNi2B2C is controlled by the static tetrahedral NiB4 
geometry. Within the model of Ref. ^3] superconductivity is more favorable for 
RNi2B2C compounds with the B-Ni-B tetrahedral angle closer to the ideal value 
of (pideai = 109.5°. For all RNi2B2C compounds (p < (pideai (for YNi2B2C (p ^ 107.3° 
jni). Then a compression along the c-axis will cause a decrease in bringing it further 
away from the (pideau and in-afe plane compression will yield an increase in (j) so that 
it will approach (pideai, consequently the signs of the uniaxial pressure derivatives of Tc 
expected from [22] are dT^/dp^ < 0, dTc/dpab > 0, opposite to that evaluated from the 
TE experiment (Table Q). 

3.2. Magnetostriction 

Magnetostriction loops taken at different temperatures, both in superconducting and 
normal state, for field applied in the ab plane, are shown in Fig. El Magnetostriction 
in superconducting state depends on number of parameters, including the extrinsic 
ones, like shape of the sample, pinning strength and its field dependence j211 ESI UEi 
|2Z|. Moreover, irreversible flux pinning induced magnetostriction causes (geometry- 
dependent) shape distortions [211 12H] that add additional requirements on mounting of 
the sample in the experimental cell for detailed study of the critical superconducting 
state via magnetostriction. Here we will just mention that in both samples the 
magnetostriction in the intermediate superconducting state is high {ALmax/ Lq > 
1 ■ 10"^), significantly higher than that seen for the polycrystalline sample in 
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both samples the features associated with peak effect j2Z| or "dip" in magnetization 
[2n] are apparently seen, with some structure for the sample B. In the normal 
state magnetostriction is rather small. In superconducting state, for magnetic fields 
appreciably lower than Hc2, the irreversible magnetostriction is much higher for the 
sample B, consistent with its lower RRR (Fig. and higher pinning in this sample. 

For field applied along c-direction the quantitative difference in the irreversible 
magnetostriction in superconducting state between two samples is even more drastic 
that for H\\ab (Fig. Very large, {ALmax/Lo ~ 6 • 10"^) effect is seen for sample B 
and is possibly associated with the peak (and/or "dip") effect. Much smaller feature in 
the irreversible magnetostriction is seen below Hc2 in the sample A. 

The magnetostriction data allow for an estimate of the Hc2 (defined here as the 
high field onset of the peak/dip effect feature). The anisotropic data for both samples 
are shown in Fig. [3 The values of Hc2 and its anisotropy {H^2/Hc2 ~ 1-05 — 1.1) are 
consistent with the previous results obtained from magnetization and magnetoresistance 
measurements [HI EUl EH 132] , being on the lower end of the reported values for the Hc2 
anisotropy. The observation that Hc2 is lower for the sample with lower RRR is broadly 
consistent with the results e.g. on Y(Nii_2.Co2:)2B2C [HSI and (Yi_2:Lua;)Ni2B2C |31j, 
however it disagrees with the generally expected increase of Hc2 with decrease of the 
mean free path (decrease of RRR) expected for conventional superconductors (see Ref. 
[HK] and refs therein). Careful study, by several techniques, of Hc2 and its anisotropy in 
RNi2B2C and in particular in YNi2B2C as a function of its mean free path is a topic of 
interest for a separate study. 

A striking feature in high field magnetostriction of the sample A is clearly seen 
quantum oscillations (Fig. IHK), that start in the still superconducting, irreversible, 
region (at approximately 0.7 Hc2 at T = 1.8 K) and continue, with growing amplitude, 
in higher fields (similarly, the oscillations are present for the sample B, but the Y-scale in 
Fig. inb does not allow to see them, see below). Quantum oscillations in magnetostriction 
is a known, albeit not so common phenomenon EZj, that requires high quality 
crystals and sensitive measurement techniques to be observed. These oscillations are 
present, with smaller amplitude, at least up to T = 10 K and are observed as well in the 
sample B. Magnetostriction data as a function of inverse magnetic field for three different 
temperatures are shown in figureslHK (sample A) andEK (sample B). At base temperature 
the amplitude of the oscillations for the sample A is ~ 4 times higher, consistent with 
lower scattering in this sample. The Fourier transformation of the data (figures |Hb, 
Eb) reveals one frequency, F = 5.04 MG. The effective mass, m*, corresponding to 
this orbit can be estimated from the slope of the ln{A/T) vs. T plot {A - amplitude 
of the oscillations in some chosen field) ^]. For samples A and B this procedure 
results in values 0.38 niQ and 0.43 mo respectively; in average: m* = 0.4 ± 0.03mo. 
Quantum oscillations in magnetization (de Haas - van Alphen effect) in YNi2B2C were 
described in number of publications (see e.g. (231331101 and references therein). The 
aforementioned oscillations in magnetostriction are consistent, both in frequency and 
effective mass values, with the strongest frequency observed by a conventional de Haas 
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- van Alphen effect measurements, the a orbit on the Fermi surface part formed by the 
17th band. 

4. Summary 

In summary, the TE measurements allow us to estimate uniaxial pressure derivatives of 
Tc which have different signs for pressure applied in the ab plane and along the c axis. 
The results call for re-evaluation of the role of NiB4 structural unit on superconductivity 
in YNi2B2C. The irreversible MS in YNi2B2C is large, complex, has features associated 
with the peak (and/or "dip") effect near and is open for further, detailed studies of 
pinning. Quantum oscillations in MS (analog of de Haas - van Alphen effect) were 
observed for H\\c. The frequency was identified as the a orbit previously seen in 
magnetization. If combined with high field magnetization measurements on the same 
crystals, the oscillatory MS will allow for estimate of evolution of the extremal orbit 
area under uniaxial stress [3^1 • 
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Table 1. Changes in specific heat and thermal expansion coefficients at Tc and 

estimates of anisotropic pressure derivatives of Tc for two samples of YNi2B2C. Sample 
A: YNi2^°B2C, solution grown, RRR « 63; sample B: YNi2^^B2C, melted zone grown, 
RRR « 16. 



Sample 






^Otab 






dTc/ dpc 


dTjdP* 




(K) 


(mJ/mol K) 


(10-« 


(10-« K-1) 


(K/kbar) 


(K/kbar) 


(K/kbar) 


A 


15.9 


444 


-7.2 


5.9 


-1.02 ■ 10-2 


0.83-10-2 


-1.21 ■ 10-2 


B 


15.1 


522 


-6.5 


6.5 


-0.74 • 10-2 


0.74 • 10-2 


-0.74 • 10-2 
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Figure 1. Temperature-dependent magnetization, resistivity and heat capacity for 
samples A and B near superconducting transition. Vertical lines mark transition as 
defined by p = 0. 
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Figure 2. Relative temperature-dependent linear and volume dilations of two 
YNi2B2C crystals. 
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Figure 3. Temperature-dependent linear and volume thermal expansion coefficients 
for two YNi2B2C crystals plotted together with literature data; symbols: Q " R-^f- IBji 
A - Ref. (see text for discussion of the literature values). 
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Figure 4. Temperature-dependent linear thermal expansion coefficients for two 
YNi2B2C crystals in the vicinity of the superconducting transition. 
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Figure 5. Longitudinal magnetostriction for two YNi2B2C crystals with magnetic 
fields applied in the a6-plane. AL is defined as L{H) — L{H = 0). 
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Figure 6. Longitudinal magnetostriction for two YNi2B2C crystals with magnetic 
fields applied along the c-axis. For sample A only data at T = 1.8 K are shown. AL 
is defined as L{H) - L{H = 0). 
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Figure 7. Temperature dependent H^z for two YNi2B2C crystals as defined from 

magnetostriction measurements. Triangles - H\\ah, circles - H\\c, filled symbols - 
sample A, open symbols - sample B. Insert shows anisotropy of Hc2- 
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Figure 8. (a) Quantum oscillations of the longitudinal magnetostriction plotted as a 
function of 1/H for sample A. (b) Fourier transform of the AL/Lq vs. 1/H curve for 
T = 1.8 K. Inset: ln{A/T) vs. T plot for amplitude A at ff ~ 133 kOe. 
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Figure 9. (a) Quantum oscillations of the longitudinal magnetostriction plotted as a 
function of 1/H for sample B. (b) Fourier transform of the AL/Lq vs. 1/H curve for 
T = 1.8 K. Inset: ln{A/T) vs. T plot for amplitude A at ~ 133 kOe. 



